Conservation laws, bright matter wave solitons and 
modulational instability of nonlinear Schrodinger equation 
with time-dependent nonlinearity * 

Shou-Fu Tian 1 ' 2 *, Li Zou 3 , Qi Ding 1 and Hong-Qing Zhang 1 . 
' 'School of Mathematical Sciences, Dalian University of Technology, Dalian 116024, 

f*^) \ People's Republic of China 

2 Department of Mathematics, University of British Columbia, Vancouver, BC, V6T 1Z2, Canada 
3 School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 1 16024, 



in 



in 



People's Republic of China 



| Abstract: In this paper, we consider a general form of nonlinear Schrodinger equation with time- 

dependent nonlinearity. Based on the linear eigenvalue problem, the complete integrability of such 
nonlinear Schrodinger equation is identified by admitting an infinite number of conservation laws. 
Using the Darboux transformation method, we obtain some explicit bright multi-soliton solutions in a 
recursive manner. The propagation characteristic of solitons and their interactions under the periodic 
plane wave background are discussed. Finally, the modulational instability of solutions is analyzed in 
the presence of small perturbation. 
PACS numbers: 02.30.Jr, 05.45.Yv, 02.30IK. 



Keywords: Exact solution, Bright matter wave soli ton, Conservation law, Modulational instability, 
Nonlinear Schrodinger equation with time-dependent nonlinearity 

^ | (Some figures in this article are in colour only in the electronic version) 

.a; 

1. Introduction 

It is known that every atom in a Bose-Einstein Condensates (BECs) moves in an effective mean 
field due to the other atoms and the mean field equation of motion governing the evolution of the 
macroscopic wave function of the Bose-Einstein condensate is the so-called time-dependent Gross- 
Pitaevskii (GP) equation [1] 
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where *F is the BEC order parameter, V ext is the external trapping potential and the coefficient g = 
4nh 2 a/m characterizes the effective interatomic interactions in the BEC through the s-wave scattering 



'Corresponding author. Tel: +86-411-84708351-8136 

E-mail address: shoufu2006@126.com, shoufu@math.ubc.ca (S. F. Tian) 



1 



length a. For a parabolic potential and time dependent scattering lengths V ext = with a{t) - 

-2zhs. t % = 2m, the above GP equation in one dimension takes the following nonlinear Schrodinger 
equation with time-dependent nonlinearity [2-5] 



Here the Feshbach-managed nonlinear coefficient a{t) can be redefined as a(f) = \a s (f)\/ciB = go exp(e?) 
(as is the Bohr radius)(6l, which is also called the time dependent scattering length. In Eq. dl.21 l, time t 
and coordinate x are measured in units 2/co ± and a±, where a ± - (h/ma> ± ) 1 ^ 2 and ciq = (h/m(Oo) ' 2 are 
linear oscillator lengths in the transverse and cigar-axis directions, respectively. a> ± and ojq are corre- 
sponding harmonic oscillator frequencies, m is the atomic mass, and the parameter e = 2\a)Q\/a) ± «: 1. 

Investigation of the behaviour of Bose-Einstein Condensates (BECs) requires solving an inho- 
mogeneous nonlinear Schrodinger equation known as the GP equation (H. Eventhough numerical 
solutions of GP equation are available QUI, construction of analytic solutions will offer more insight 
into the domain of BECs opening the doors for developing concrete applications of BECs in future. 
As we well know, it is significantly important in mathematical physics to search for exact solutions to 
equation dl.2| ). Exact solutions play a vital role in understanding various qualitative and quantitative 
features of nonlinear phenomena. It is well known that searching for soliton solutions of the nonlin- 
ear evolution equations is one of the most important topics in soliton theory. Darboux transformation 
(DT) (9J [10l has been proven to be one of the most fruitful algorithmic procedures to obtain exact 
solutions of the nonlinear evolution equations. 

The main aim of the present paper is to construct some infinite number of conservation laws, ex- 
plicit bright multi-soliton solutions by using DT method, and investigate the modulational instability 
of solutions of a general form of nonlinear Schrodinger equation with time-dependent nonlinearity 
(1 1 -2b - In this paper, on the basis of the Lax pair associated with Eq. (11.2l ). we will derive an infi- 
nite number of conservation laws to identify its complete integrability. Furthermore, we will apply 
the Darboux transformation method to this integrable model and give the general procedure to re- 
cursively generate the bright jty-soliton solutions from an initial trivial solution. Moreover, we will 
discuss the propagation characteristic and interactions of solitons under periodic plane wave back- 
ground and analyze the linear stability of the nonlinear plane waves. 

The paper will be organized as follows: In Section 2, based on the linear eigenvalue problem 
associated with Eq. (ll.2b and an auxiliary ID NLSE (12.21) , we obtain a series of conservation laws. 
Then the integrability is identified by admitting an infinite number of conservation laws. In Section 
3, we will apply the Darboux transformation method to this integrable model and give the general 
procedure to recursively generate the bright A^-soliton solutions from an initial trivial solution. The 
propagation characteristic of solitons and their interactions under the periodic plane wave background 
are discussed in Section 4. In Section 5, we analyze the modulational instability of the nonlinear plane 
waves. Finally, some conclusions and discussions are provided. 
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2. Lax pair and infinite conservation laws 



In this section, using the linear eigenvalue problem associated with ID NLSE (11.21 ) and an auxiliary 
ID NLSE 02.21 ), we construct an infinite number of conservation laws for the ID NLSE (I1.2I ). Then 
the integrability is identified by admitting an infinite number of conservation laws. 
Under the following transformation 



ifr(x,t) = u(x,t)e 't f 



(2.1) 



to Eq. (11.2l) . we can obtain the following new ID NLSE 

d d 2 d 

i—u(x, t) -\ ~u{x, t) - iex—u(x, t) + 2a(t)\u(x, t)\ 2 u(x, t)e~ et - ieu(x, t) - 0. (2.2) 

ot dx 1 ox 

By virtue of the Ablowitz-Kaup-Newell-Segur scheme |[TTl[T2l . the Lax pair associated with equation 
( 12.21 ) can be derived as 

<D. t = U®, <D f - V<D, (2.3) 

where $ = (fa, fa) T (the superscript T denotes the vector transpose) is the vector eigenf unction (2], 
and the matrices U and V have the following forms 
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Q = 





U = 


AJ + P, 






V = 


2iA 2 J + eAxJP + 2ieP + Q, 
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K -^gou(x,i)* 






igo\u(x, t)\ 2 - ^exu(x, t) + i y/gou(x, t) x 


k ^exu(x, t)* + i yffftu(x, t)* x -ig Q \u(x,t)\ 2 



(2.4a) 
(2.4b) 



where go - a(t)exp(-et) is an arbitrary function. Hereafter the asterisk stands for the complex 
conjugate. From the compatibility condition U t - V x + [U, V] = 0, one can derive Eq. (l2.2b and 
Eq. dl.2l ) in the case of tf/(x, t) = u(x, t)e~ l ~~2, respectively. 

Using Lax pair (12.31) and Refs. lfT3"l[T4"1 . we can further derive an infinite number of conservation 
laws of Eqs. (ll.2b and (12.2b - By introducing the quantity f(x, t) = A/gou(x, t)^r, the linear equations 
(12.31 ) can be transformed into the following Ricatti equation 



Ui^JC^ t^) XX p U\JCy t/X t*/ \ r\ * /*/■ \'2. n * p i~\ 

f(x, t)x + i -fix, t) + a(x, t) f(x, t) — 2igof(x, t) — 2igQii(x, t)u(x, t) + j3(x, t)f(x, t) = 0, 



u(x, t) 



u(x, t) 



(2.5) 

with a(x, t) - 2ie - 2i + eAx - ex + iA and yS(x, t) = eA - e + eA x + 2ieA - eAx - 4iA . Substituting 
f(x, t) = Z^li Jtjtf * nt0 ^i catt i equation (12.51 ) and equating the like powers of f5(x, t) to zero, we have 
a recursion formula 



f\ - 2ig Q uu*, f n+l = 2ig V fif n _i - f nx - i—f n - a—f n = 0, (n = 2,3, • • • ), 

4— 1 u u 

1=1 



(2.6) 
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where /„ (n = 1,2, • • • ) are the functions to be determined. By virtue of the compatibility condition 
(In cpi) rf = (In (f>i) fr , we obtain the following conservation form 

i- Pk (x,t) + — 3*(jc,0 = 0, (2.7) 
at ox 

where pt(x, t) and 3fc(x, (k - 1, 2, • • • ) are called conserved densities and conserved fluxes, respec- 
tively. The first three significant physical conservation laws are presented as 

Pi(x, t) = 2/goM 2 , Pi(x, t) = -2ig (\ + a)u x u* - 2ig uu* + 2g u xx u*, 

P3,(x, t) = -8/,£>q|m| 4 - 2gou xx u* + 2/(1 + a)u xx u* + 2iego(A - \)u x u* + 2/(1 + 3a)gou x u* 

[2/go^:t: M * _ 2ia{\ + a)gou x u* + 2(1 + 2a)gou x u xx u*] , 

3i(x, t) = /goM 2 - 2g Vio«x"* + 2i(eAx + lie - ex)g Vgo««\ 

2 i ^ U % U ^ , Liy.lt . 

3 2 (x, t) = igo\u\ + 2g yfgou x u x + 2ig ygo ■ + 2(1 + a)g V#o— 2ig yfgoieAx + lie 

— ex) [uu* x + iu xx u* + (1 + a)u x u*] , 

2 / u x \ 

33(x, t) - igo\u\ + ^Jgo \eAx + lie - ex + / — j 

X [-8/£q|m| 4 - lgQU xx u* x + 2/(1 + a)u xx u* + 2ie{A - l)g u x u* + 2/(1 + 3a)go«x M ^] 

— V#o + _ ex + i—\ \lig u xx u* - 2ia(\ + a)gou x u* + 2(1 + 2a)gQU x u xx u*] . 

LI \ 11/ 

The conserved quantities J! = 2igo J_™ \u\ 2 dx, J%? = -2igo f J^O + a)u x u* + uu* x + iu xx u*dx and 
,J4f = 2igo -Ag^\u\^+iu xx u x +^{\+a)u xx u* +e{A-l)u x u* +{\+3a)u x u* x -j t [u xx u* - a(l + a)u x u* 
-/(l + 2a)u x u xx u*] dx represent the energy, momentum and Hamiltonian, respectively. 



3. Darboux transformation and bright soliton solutions 

In this section, the Darboux transformation method is applied to this integrable model and the general 
procedure is presented to recursively generate the bright A^-soliton solutions from an initial trivial 
solution. Based on Lax pair (12.31) . the use of the Darboux transformation to construct the bright soliton 
solution of nonlinear partial differential equations is an optimum choice among various methods in 
soliton theory [15-26]. Owing to its purely algebraic algorithm, with the symbolic computation, 
the analytical N-soliton solution can be generated through successive application of the Darboux 
transformation [31,32]. 



3.1 Darboux transformation 

The Darboux transformation is actually a gauge transformation 0[1] = T$> of the spectral problem 
(12.31 ) by considering the following gauge transformation 

<D[1] - (AI - S)<D with 5 = HAH' 1 , A = diag(^i,^), (3.1) 

where H is a nonsingular matrix. It is required that <1>[1] solves the same spectral problems ( 12.31 ) 

0[l] v = £/[l]<D[l] and <D[l] f = V[1]0[1], (3.2) 
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with U[l] = AJ + P[l], V[l] = 2iA 2 J + eAxJP[l] + 2ieP[l] + Q[l] and 

i 

P[l] = 





' 1 ^ 


, / = 
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k - 1 > 



em 



V^o" 

{-^u[\r o 

^ol«[l]| 2 - yfgaexu[l] + i Vgo«[l];c 
y ^exu[l]* + i*JgEu[l]* x -ig Q \u[l]\ 2 



(3-3) 



where u[l] = i//[l]e' 4 + 2 . Substituting Eq. (13.11 ) into Eqs. (I3.2I ). we have the following relationship 

P[l] =P + [/,5], 5, + [S,JS +P] = 0. (3.4) 

Now we discuss a concrete transformation. It is easy to verify that if O = (cj)i , fa) 7 is an eigenfunction 
of Eqs. (12.31 ) with eigenvalue A = X\, then {fa, -fa,) T is also an eigenfunction of Eqs. (12.31 ) with 



(3.5) 



eigenvalue A = A* x . Thus we take the matrix H in the form 

fa (j>* 



1 4> 2 j 

Therefore, by means of Eqs. (13.41) and (13.5b . the onceiterated new potential of Eqs. (ll.2b and (12.21 ) are 
given by 



tA[l] = u[Y[e 
u[l] - u + 2 



(Ai-A^fafa 



(3.6) 



fa$\ + fa§2 

It is straightforward to verify that the Darboux transformations (13.11 ) and (13.61 ) can simultaneously 
preserve the form of linear eigenvalue problem (12.31) . 

According to the above computation, we have the following propositions. 



Proposition 3.1. When u(x, t) is given, let {fa, fa) T be the solution of Eqs. ft23i with A = A^, then by 
using of the Darboux matrix T (13.11) and Darboux transformation (13 -6b . we have 



U[l] = (T x + TU)T-\ with T = AI-S. 
We can obtain the same proposition about the auxiliary spectral problem. 



(3.7) 



Proposition 3.2. When u(x, t) is given, let (fa-, fa) T be the solution of Eqs. (12.31 ) with A = Au, then by 
using of the Darboux matrix T (13.11 ) and Darboux transformation ( 13.61 ), we have 

V[l] - (T, + TV)T' 1 with T = AI - S. (3.8) 

Proof. Substituting Eqs. (13.11 ). (13.21 ) and (13.61 ) into Eq. (13.81 ) by a direct calculation, we can obtain 
the conclusion. The proof is completed. □ 
Propositions 3.1 and 3.2 show that the transformation (13.11 ) and ( 13.61 ) change the Lax pair (12.31 ) 
into another Lax pair of the type (13.21 ) with U[l] and V[l] having the same form as U and V, respec- 
tively. Therefore both of the Lax pairs lead to the same equation ( 12.21 ), so Eqs. (13.11 ) and ( 13.61 ) are the 
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Darboux transformation of Eq. d2.2l ). From propositions 3. 1 and 3.2, we have following theorem. 



Theorem 3.3. Assuming that if/, u, A\\\, c|>[l, Ai\ 2 ) T , (<t>[2, A 2 ]i, cb[2, A 2 ] 2 f , ■■■ , (§[N, A N ] U §[N, A N ] 2 ) T 
be the solution of the ID NLSE (12.21 ) and N linearly independent solutions of the linear eigenvalue 
problem (12.3b . respectively, and after iterating the Darboux transformation (13.11 ) and (13.61 ) ./V times 
analogous to the above procedure, we can further obtain the Nth-iterated potential transformation as 

2 

if/[N] = u[N]e- iM r-?, 

u[N] = u + 2 > ■ , (3.9) 

with 

§[i + l,A i+1 ]j = (A M - A*)dp[i, A i+1 ]j - ~^(h ~ A*)$[i, /I;]/, 
sdi = dp[i,Ai]\§[i,Ai + i]i + <$>[i, AiY 2 $[i, A i+1 ] 2 , 

3§i = #i, A i+1 ]mai + l]l + W,Ai +l ] 2 ^[i,Ai + iT 2 , d= 1,2,--- ,N-\,j=\,2), (3.10) 

where {<§>\k, A^\\ , <$>[k, Ak\ 2 ) (k = 1, 2, • • • ,N) is the eigenfunction of Eqs. (12.3b with the eigenvalue 
A = Ak and potential u[k - 1]. 

Proof. (By induction) The case N = I follows from Darboux transformation (13.11 ) and (13.6I ). Now 
assume that equation (13.91) holds for N = n. Then 

• eA 2 el 

if/[n] = u[n]e 4 2 , 

u \n \ = u + 2 > , (3.1 1) 

^ <\>[i, Aih<$>[i, A^ + AihtyV, A { ]* 

where 

())[/ + 1, j = (Am - A { A i+ i]j - —(Aj - A { Ai] j, 
sii = AiY^U, A i+l h + AtY 2 fy[i,A i+l ] 2 , 

&i = <\>[i,A i+l ] l( \>[i,A i+ n\+W^i + ih<\>[i,Ai + iY 2 , (i = 1,2,- •• ,n - l,j = 1, 2). (3.12) 



In case of N = n + I, using Darboux transformation (13.61 ). one obtains 

(A, 1+l - A* n+l )$[n + 1, A.„+i]icj)[?2 + l,l„ +l Y 2 



u[n + 1] =u[n] + 2 



<\>[n + l,A,„ + i]ict>Di + l,A,„ + i]* + c()[n + 1, X.„ +1 ] 2 cb[n + l,l n +\Y 7 



n+l 



(Aj - A*)d>[i, Ai]i&>[i, AiY, 
^ (()[;, /!,] i Ai]\ + (\>[i,Ai] 2 §[i,AiY 2 

Based on the first equation of Darboux transformation (13.6I ). it is easy to obtain 

+ 1] = m[« + l]e" ! V-!r. (3.14) 

From Eqs.(f3ll3]) and (l3~T4l ). the conclusions of Eqs.(f3S and (l3~T0l) are hold for N = n + 1. This 
completes the proof. □ 
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3.2 Bright matter wave soliton solution 

Bright soliton here is defined as the emergence of a positive pulse. Bright soliton solutions have been 
investigated by original Ref s. ifTTl |27l |28l . In the following, the Darboux transformation is applied to 
construct the explicit bright soliton solutions of Eq. dl .2b . Using the trivial solution - (u=0), we 
solve the linear equations ( 12.31 ) with A = Ai = +iv\) and obtain the eigenf unction 

^ = J(Pi+i^) t fa = e -fa+Hh) t (3.15) 

with 



Pi 



■ fiix- 2/zivi^i, #i = vix + fal -vf)t + t\u (3.16) 



where "E,\ and r|i are arbitrary constants. 

The bright one-soliton solution for Eq. d 1 .2b are obtained as 



: or er 
2 



\jj(x, t) - u(x, t)e 1 4 
u(x, t) = i ' v i e 1 sechp i . (3.17) 



by substitution of the above results into formula (13.6b . It implies that the imaginary part v\ of the 
eigenvalue A\ determines the amplitude of the solitons \fr(x, t) and u(x, t), while the velocity of solitons 
are related to both real and imaginary parts of the eigenvalue A\. 
From above, we have the following propositions. 

Proposition 3.4. Assuming that (c|)[l, A\\\, §W,A\\2) T and (cj)[2, ^h, <\>[2, A2]2) T be two linearly in- 
dependent solutions of the linear eigenvalue problem ( 12.31 ) corresponding to two different eigenvalues 



A\ - \{p\ + iv\) and A2 - \(pLi + iv%), respectively, we obtain the bright two-soliton solutions of 
EqAl.2j and EqA2.2j from formula ( 13.91 ) with N - 2 

i/f[2] = u[2]e' i ^, 

^e^coshpi + ^ 2 e i,?1 cosh^2 + % (e i,?1 sinhp2 - e^ 2 sinhpi) 

i (3 18) 

^cosh(pi + P2) + "^cosh(pi - P2) - 2viV2Cosh(??2 - #1) 

with 

pi = p t x - iMVi&.i, = v t x + (ju? - vfj t + r\i, 

% = v 2 [(mi -m) 2 ~ v? + vl] , %, = vi [Oui -p 2 ? + y\ - vl] , 

1 , 1 2 

^3 = -2iviv 2 (pi - Mi), Vl = rO"l ~M2) + 2^ Vl ~ v ^> ' 

r 2 = ^(Mi-M2) 2 + j(vi+v 2 ) 2 , (i = 1,2), (3.19) 
where %i, T],- are all complex constants and the parameters pi, Vj + 0(i = 1,2). 

Proof. It is straightforward to prove this proposition by using the initial values (13.151 ). (13 - 1 6b . and 
Theorem 3.3 for u = 0, N = 2. □ 
Proposition 3.5. Assuming that (4>[1, X{\\, $[1, A\\ 2 ) T , ($[2, ^2] 1 , <\>[2, A 2 h) T , • • • , ^wh) 7 
be TV linearly independent solutions of the linear eigenvalue problem (12.31 ) corresponding to N differ- 
ent eigenvalues A\ = j(Mi + ivi), ^2 = + ivi), ■ ■ ■ ,An = j(pn + ivn)> respectively, we obtain the 
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bright N-soliton solutions of Eqs. ( fL2| ) and ( 12.21 ) from formula ( 13.91 ) with N 



iff[N] - MlAfle-'V-T, 



i[iV] = 2j] 



1=1 



cj)[j, /li]i<|)[i, A,]! + <j)[i, ^] 2 4)[/, Ai\* ' 



with 



pi = mx - IpiiV&i, &i - v t x + (ju? - v?) / + T];, (i = 1, 2, • • • ,iV), 



(3.20) 



(3.21) 



where n,- are a/Z rea/ constants and the parameters /Uj, v; ^ 0(/ = 1 , 2, • • • , iV). (c()[^, i , §[k, Ak]t) T 



satisfies Eq A3.l0i and is the eigenfunction of Eqs. ( 12.31 ) with the eigenvalue A - and potential u[k] 
(k = 1,2,- •• ,N). 

Proof. It is straightforwai - d to prove this proposition by using the initial values (13.151 ). (13.161 ). and 
Theorem 3.3 for u = 0. □ 
The graphs of the bright one-soliton periodic wave solution (13.171) and two-soliton periodic wave 
solution ( 13.181 ) are plotted in Fig. 1 and Fig. 2, respectively. 





(c) 

Fig. 1. (Color online) A symmetric bright one-soliton solution \if/(x, f)| 2 of Eq. d 1 .2b with parameters: yUj = 1, 
Vi = 0. 1, ^i = 1 + i, T]i = 3 and e = 0.0 1 . This figure shows that the symmetric bright one-soliton periodic wave 
is spatially periodic in two directions, but it need not be periodic in either the x or t directions, (a) Perspective 
view for the real part of wave, (b) Perspective view of the imaginary part wave, (c) Wave propagation pattern 
of the real part wave along the x axis, (d) Wave propagation pattern of the imaginary part wave along the x 
axis. 
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(c) (d) 
Fig. 2. (Color online) A symmetric bright two-soliton solution \if/(x, t)\ 2 of Eq. d 1 .2b with parameters: fi\ = 0.1, 
v\ = 0.1, |i = i, r\i = 0.1/, jU2 = 1, vi = 1, 12 = i, i>]2 = 0.1;, and e = 0.01. This figure shows that the symmetric 
bright two-soliton periodic wave is spatially periodic in two directions, but it is periodic in either the x or t 
directions, (a) Perspective view for the real part of wave, (b) Perspective view of the imaginary part wave, (c) 
Wave propagation pattern of the real part wave along the t axis, (d) Wave propagation pattern of the imaginary 
part wave along the t axis. 



4. Bright soliton interactions on the periodic background 



In this section, we investigate some bright soliton interactions on the periodic background by using Eqs.(|37 
and (13.9b . The simple exact solution of Eq. dl .2t is the plane wave 



i//(x, f)o = u ( x , t)a e ' 4 



u(x, f)o = ye 



Hkx+mt) 



(4.1) 



where y and k are real constants, and the frequency to solves the nonlinear dispersion relation 

03 = -k +ek + 2y - z'e + — . 

7 4 



(4.2) 



Considering solution ( 14.1b as the initial seed solution of Eq.( fT72b , we obtain the linear equations ( 12. 3t and have 
the eigenfunction corresponding to the eigenvalue Ai in the form 



cMMih =d\e Bi +d 2 e s ' 1 , o>[l, Ai] 2 = d 3 e~®< + d 2 e~® 2 



(4.3) 



with 



^ = -i(kx + rof) + T(x + At), 



$ 2 = -i(kx + cof) - T(x + At), 



ik 



T - A] 



u ' + 5)-^*=^b(l +r - J ' 



A = ^/li - ek - 2ieA\ - i— — — + iy - i— + {lie - k) 



/li + - 



(4.4) 
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where d\ and di are two arbitrary complex constants. 

Substituting expressions (14.3b into formula d3.9t , the one-soliton and multi-soliton solutions on the pe- 
riodic background can be obtained with the iterative algorithm of the Darboux transformation. The bright 
one-soliton solutions are plotted in Figure 3, which shows four kinds of soliton profile structures with different 
wave numbers. Although solitons in Figures 3(c) and 3(d) hold larger wave numbers than those of Figures 3(a) 
and 3(b), they both can propagate stably for long distances by the results of numerical simulations of nonlinear 
pulses. 



Proposition 4.1. Assuming that (cj>[l, cp[l, be the solution of the linear eigenvalue problem (12. 3\ 

with eigenvalues A\, we obtain the bright two-soliton solutions of Eqs. (II. 2\ and d2.2l ) corresponding to eigen- 
values A2 from formula d3.9t with N — 2 



ij,[2}=u[2]e-'---, 

2 

u[2] = u(x,y) Q + 2 2^ 



(^ - Ap#Ui]i<t>[U,] 



The solution of the linear eigenvalue problem (12.3b with eigenvalues A2 is 



W2M 
cj>[2,A 2 ] 2 



(A, + AIMUAM^hAM 2 - A'MUA^UAi]; - A^hA^MhMt 



(j)[l,Al]l(t)[l,^l]* + #l,/ll]2Ct)[l,/ll]* 
■^)|#l,Ai]i| 2 #l,Ai] 2 + ^|<|)[l,Ai]i| 2 <|)[l,yli]J-Al<|>[l,yli]| 



WhAnMhAAi + whMhWhM]; 



(4.5) 



(4.6) 



The bright two-soliton solutions are plotted in Figures 4. Figure 4 depicts the bright two-soliton periodic 
wave solutions in the three-dimensional space. From the interaction process in two sets of figures, it can be 
clearly seen that the interacting solitons like particles cross each other unaffectedly only by a phase shift, and 
their respective amplitudes and velocities are the same as those before collision. Further theoretical analysis 
for soliton solutions shows that the sign and value of real part of the eigenvalue determines the propagation 
direction of soltion and the amplitude, respectively. 
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Fig. 3. (Color online) Four kinds of bright one-soliton solutions \if/(x, t)\ 2 of Eq. jl.2b with different wave 
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numbers k. The related physical quantities are y — 2,d\ — \,di — 3,Ai — 5 - 2;, go = 3,e = 0.001 and wave 
numbers (a) k = 1, (b) k = 10, (c) & = 50, (d) k = 100. 
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Fig. 4. (Color online) A asymmetric bright two-soliton solution \iff(x, t)\ 2 of Eq. dl.2b with parameters: y = 
2,d\ = I, do = 3,/li = 10 + 2;, /I2 = 10 - 2i, go = l,e = 0.001 and wave numbers k = 1. This figure shows 
that the asymmetric bright two-soliton periodic wave is spatially periodic in two directions, but it need not be 
periodic in either the x or t directions, (a) Perspective view of the wave, (b) Wave propagation pattern of the 
wave along the x axis, (c) Wave propagation pattern of wave along the t axis. 



5. Linear stability analysis 

In this section, we analyze the modulational instability of the nonlinear plane waves. The phenomenon of mod- 
ulational instability of nonlinear waves has a significant importance in the theory and experiment of nonlinear 
waves. The reason is that this phenomenon is quite relevant applications in many different branches of physics 
such as in condensate physics, plasma physics, hydrodynamics and nonlinear optical fibers l29l . Of particular 
interest for studies of stability of the solution of the NLS equation is that the envelope of a monochromatic 
plane wave propagates in nonlinear medium in the presence of microwave noise pulses. 

As we well-known that the exact plane wave solution for the focusing NLS equation is to be linear unstable 
l30l . According to the method of investigating the linear stability of the NLS equation, we perform a linear 
stability analysis of the nonlinear plane waves solution for Eq. Jl .2t . The solution ( 14. U is the plane wave with 
constant amplitude for Eq.( ll.2l ). We investigate the solution with the perturbation in the form 

\j/{x, t) = u(x, f)e~'T~~?, 

u(x, t) — uq(x, f) [1 + eu(x, t)] . (5.1) 
By substitution of ( 15. U into Eq. d2.2t , the linearized disturbance equation at o(e) becomes 

iu t + u xx + i(2k - xe)u x + 2y 2 a(t)e e '(u* + 2ft) + e(x - 1) ( -x + - + k) ft = 0. (5.2) 
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By virtue of the linearity of Eq. d5.2b , the solution can be expressed as the following linear combination 



u(x, t) = ^ + e mx - @,) + S£le 



-W{x-Wi) 



(5.3) 



We characterize this solution by the real disturbance wave number and the complex phase velocity 0. Sub- 
stituting ( 15.31 into d5 .2b and collecting resonant terms, we have two linear homogeneous equations 



0@-0 2 - 0(2k - xe) + 4y 2 a(t)e a + e(x - 1) f |jc + | + kj + 2y 2 a(t)e e 'Jtf- = 0, 



2y 2 a(t)e"^ + + 



-0® -0 2 + 0(2k - xe) + 4y 2 a(t)e a + e(x - 1) f |* + | + kj 



(5.4) 



By using the determinant of the matrix of coefficients of linear Eqs. ( 15.4b . the dispersion relation for linearized 
disturbance can be expressed as 



© =2k -xe± 
-2k — xe ± 



4y 2 a(t)e et + 2 



-(H 



- 4y 4 a(t) 2 e 2el 



i + kj" \4y 2 g a e 2 « + 02 - e (1 + kjj - 4y 4 g 2 e^ (with a(t) = g e et ), (5.5) 



where < e « 1 and y, k are real constants. 

By virtue of the above relation ( 15.51 ), if satisfies one of the inequalities as follows 



the frequency is real at any value of the wavenumber 0, whereas becomes complex. 



(5.6) 
(5.7) 



Proposition 5.1 (i) For Eg. d5.61 l, if 2 > e(| + kj and k > i > i-|(0<e<sc I), the frequency is real at 
any value of the time variable t. Otherwise, there is the instability region of the modulational waves, and the 
disturbance will grow with time exponentially, (ii) For Eq. (5. 71 ), ifO < 2 < ~e(f + kj and k < — |(0<e<*: 1), 
the frequency is real at any value of the time variable t. Otherwise, there is the instability region of the 
modulational waves, and the disturbance will grow with time exponentially. 

Compared with the analytical method, in the near future, we intend to study the modulational instability 
numerically for the ID NLSE dl.2t since the linear stability analysis actually shows that this instability is 
possble. 



6. Conclusions and discussions 

In this paper, we investigate a general form of nonlinear Schrodinger equation with time-dependent nonlin- 
earity, which describes conservation laws, bright matter wave solitons and modulational instability in Bose- 
Einstein condensates with the time-dependent interatomic interaction in an expulsive trapping potential. By 
virtue of the Lax pair, we have further confirmed the complete integrability of such a model by deriving an 
infinite number of conservation laws. We have adopted the Darboux transformation method to construct the 
multi-soliton solutions via algebraic iterative algorithm. We have discussed four kinds of bright one-soliton 
solutions \if/(x, t)\ 2 with different wave numbers k, and two-periodic bright solitons of Eq. dl.2l ) under the pe- 
riodic background. Finally, we have analyzed the linear stability of nonlinear plane waves in the presence of 
small perturbation. In addition, it can be also applied to other spectral problems in the field of nonlinear partial 
differential equation appear in mathematical physics. 
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